Quantum phases of a frustrated four-leg spin tube 
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We study the ground state phase diagram of a frustrated spin-1/2 four-leg tube. Using a variety 
of complementary techniques, namely density matrix renormalization group, exact diagonalization, 
Schwinger boson mean field theory, quantum Monte-Carlo and series expansion, we explore the 
parameter space of this model in the regime of all-antiferromagnetic exchange. In contrast to 
unfrustrated four- leg tubes we uncover a rich phase diagram. Apart from the Luttinger liquid 
fixed point in the limit of decoupled legs, this comprises several gapped ground states, namely a 
plaquette, an incommensurate, and an antiferromagnetic quasi spin-2 chain phase. The transitions 
between these phases are analyzed in terms of total energy and static structure factor calculations 
and are found to be of (weak) first order. Despite the absence of long range order in the quantum 
case, remarkable similarities to the classical phase diagram are uncovered, with the exception of the 
icommensurate regime, which is strongly renormalized by quantum fluctuations. In the limit of large 
leg exchange the tube exhibits a deconfinement cross-over from gapped magnon like excitations to 
spinous. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.10.Dg, 75.10.Kt, 



I. INTRODUCTION 

Quasi one-dimensional spin systems, comprising chain, 
ladder and more involved magnetic structures are an ac- 
tive field of research thriving on a constant feedback 
between material synthesis, experimental investigations 
and theoretical predictionai"— . 

Magnetic frustration is a key issue in this field, which 
has experienced an upsurge of interest, starting with 
the discovery of J1-J2 chain materials, like CuGeOg^, 
followed by the investigation of spin tube compounds 
with an odd number N of sites per unit cell, such as 
[(CuCl2tachH)3Cl]Cl2^ and CsCrFzf^ with iV = 3, and 
^a.2^?,0'^ with iV = 9. Spin tubes with an odd num- 
ber of legs and only nearest neighbor antiferromagnetic 
(AFM) exchange are geometrically frustrated. Because of 
the Lieb-Schultz-Mattis theorem the ground state of such 
systems is either gapless and non-degenerate, or gapped 
with a broken translational invariance. Indeed, for spin- 
1/2 tubes with A'^ = 3 a spin gap was found in case of 
identical couplings on the triangular rungs, with a tran- 
sition into a gapless and translationally invariant phase 
at already weakly non-equivalent couplings.—"— N = 3 
spin-1/2 tubes with isosceles triangle basis also show a 
1 /3 magnetization-plateau 

Recently, Cu2Cl4-D8C4S02 has been established as 
a new spin-1/2 tube with an even number of legs^^, 
namely A^ = 4. Tubes with A^ = 4 and only near- 
est neighbor AFM exchange are not frustrated. How- 
ever, substantial next-nearest neighbor AFM exchange, 
diagonally coupling adjacent legs, has been claimed for 
CU2CI4 -0804802, rendering also this ladder system frus- 
trated. Inelastic neutron scatterin g^^'^^ has revealed 
a strongly one-dimensional (ID) elementary excitation, 
which is gapped and slightly incommensurate. The for- 
mer is consistent with Haldane's conjectureii for ID spin 
systems with an even number of spin-1/2 moments per 
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Figure 1: (color online) (a) Frustrated four-spin tube. Solid 
circles represent spin-1/2 moments. Plaquettes (bold black 
lines) are coupled by nearest (Ji) and next nearest (J2) an- 
tiferromagnetic exchange, blue and red lines, respectively. 
On- plaquette coupling is Jo. (b) Frustrated four-spin tube 
unwrapped, displaying structure of an anisotropic triangular 
lattice on a torus. 



unit cell. The latter is consistent with a frustrated ex- 
change. Magnetic fields have been shown to stabilize the 
incommensurate spin correlations 

Motivated by this, a geometrically frustrated and sim- 
plified four-spin tube (FFST) model has been introduced 
in Ref . M 
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with a lattice structure and exchange couplings Jim as 
shown in Fig. [T] Spin-1/2 moments are located on the 
solid circles and all couplings Jo, 1,2 are antiferromagnetic 
(AFM). The FFST lattice is identical to an anisotropic 
triangular lattice on a torus with four site circumference. 

For Ji^2 ^ Jo, the quantum properties of the FFST 
can be understood in terms of weakly coupled four-spin 
plaquettes, which allows for series expansion in terms of 
Ji^2- In Ref. [1^ such a series expansion has been carried 
out in detail regarding the one- an two-particle excita- 
tions in this restricted parameter regime. However, an 
understanding of the quantum phases of the FFST on a 
larger scale is still missing. 

Therefore, in this paper we will present combined re- 
sults from a large variety of complementary methods, 
namely, density-matrix renormalization group (DMRG), 
exact diagonalization (ED), scries expansion (SE), 
Schwinger bosons mean field theory (SBMFT), and quan- 
tum Monte-Carlo (QMC) in order to explore the param- 
eter space of the FFST. We set Jo = 1, except where 
explicitly indicated, and denote by L the tube length. 

The structure of the paper is as follows. In section HIl 
we briefly summarize the phase diagram of the classical 
FFST. In section Hill we consider the quantum case focus- 
ing the discussion onto the strong and intermediate on- 
plaquette exchange in subsection IIII Al and on the limit 
of very large leg exchange in subsection IIII Bl Section Hvl 
summarizes our picture of the FFST. For completeness 
we briefly summarize some of the methods used and refer 
to important references for them in appendix VK\ 



II. CLASSICAL PHASE DIAGRAM 

While we are studying a quantum model in one di- 
mension which does not allow for breaking of a con- 
tinuous symmetry at zero temperature, it is neverthe- 
less instructive to compare the quantum case consid- 
ered in the following sections with a classical magnetic 
phase diagram of the FFST. From Ref. [3 it is known 
that there are four ordered regimes in which S(r/) = 
S'(cos(Q • r;), sin(Q • r;), 0), where r/ is a lattice site: 

1- J2 < (1 + 2Ji)/(2(Ji + 1)) and J2 ^ Jf. commen- 
surate Q = (tt, tt) AFM 

2. J2 ^ (1 - 2Ji)/(2(Ji - 1)), with Ji < 1, and J2 ^ 
Ji: commensurate Q = (0,7r) AFM 

3. J2 ^ (2Ji - l)/(2(Ji - 1)), with Ji > 1: commen- 
surate Q = (tt, 0) AFM 

4. Two degenerate incommensurate spirals with 
Q(Ji, J2) = ±(2 arctan(Q!), 7r/2), and a — (Ji + 
\J Ji + J2) / J2 in the remaining region. 

These are shown in Fig. [21 All classical transitions are of 
first order. 
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Figure 2: (color online) Four classical phases of the FFST 
with commensurate pitches Q = (11,11), (7r,0), (0, vr) and in- 
commensurate regime Q(Ji, J2) as in text. 



III. QUANTUM PHASES AND CORRELATION 
FUNCTIONS 

In the following we gather information from various 
complementary methods to develop a quantum version 
of the phase diagram of the FFST. The dicussion is split 
into two subsections. The first focuses on the strong 
and intermediate on-plaquette exchange and comprises 
an analysis of the ground state energy using DMRG, SE, 
SBMFT, and ED, followed by an evaluation of the phase 
diagram from SBMFT, and finally a DMRG study of 
correlation funtions and structure factors. In the second 
subsection we analyze the spin excitations in the limit 
Jo ^ using QMC. 

To begin, we note that in the quantum case and at 
the points Ji(2) — ^ 00, J2(i) = the FFST is in a 
Luttinger liquid (LLQ) state. Staying on either of the 
two axes ( Ji(2) 7^ 00, J2(i) = 0) , the system is unfrus- 
trated, the inter-leg coupling is relevant, and the FFST 
opens a spin gap. This gapped phase is adiabatically con- 
nected to that of unfrustrated weakly coupled plaquettes 
( Ji(2) 1, J^( n = 0) which has been studied extensively 
in Refs. [l9ll20l The frustrated weakly coupled plaquette 
regime shows no transition between a (tt, tt) and (0, tt) 
phase, rendering the diagonal line in the lower left cor- 
ner of Fig. 121 a classical-only effect. 



A. Strong and intermediate on-plaquette coupling 

1. Ground state energy 

A natural question arising is, how far the weakly cou- 
pled plaquette phase extends away from the Ji(2) axes 
lines and if its break down is of first or second order. We 
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Figure 3; (color online) Ground state energy per site eo for 
Jo — 1. Panel (a) eo vs. Ji at J2 = 0. Panel (b) eo vs. 
J2 — Jl. Green solid: plaquette SE at 0(7). Dotted blue 
(red): ED with L=6 and PBC (SBMFT with L=400 and 
PBC). Large crosses: DMRG with L=20, ni=300, and ORG. 
Inset in (b): first order transition line in the Ji,2-plane. 



check this in two ways, considering the ground state en- 
ergy eo versus Ji_2 and the static structure factor. The 
results for eo are summarized in Fig. [3l It depicts the re- 
sults from different techniques (appendix IX)) . along two 
paths in parameter space. Panel (a) is along the Ji-axis, 
while panel (b) diagonal path Ji = J2 of maximum frus- 
tration. 

Along the Ji-axis, panel (a) the energy is a smooth 
function. All methods are in satisfactory agreement up 
to Jl « 0.7. At this point the bare SE shown, which 
has been obtained up to 0(7) (appendix I A loses con- 
vergence, while the other techniques continue to agree 
throughout the range shown. We note that finite size 
effects on the DMRG and ED are expected to be small 
since the system is gapped. 

Along the line of maximum frustration, Fig. [3] (b), the 
energy as obtained from DMRG and ED shows an obvi- 
ous discontinuity in its first derivative at Ji « 1. This 
signals a first order quantum phase transition. Remark- 
ably this point is rather close to the classical tricritical 
point, separating (7r,7r). (0,7r) and spiral classical phases 
of Fig. [21 By construction SE based on a single unper- 
turbed starting state is unable to detect this transition, 
which is consistent with Fig. [3] (b), where the SE agrees 
perfectly with DMRG and ED exactly up to the kink in 
eo. Finally SBMFT is very close to DMRG and ED in 
this panel beyond the transition, however it undcrcsti- 
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Figure 4: (color online) SBMFT phase diagram. Pitch vectors 
label short range spin correlations. Grey regions correspond 
to unphysical 'decoupled chain' phases 



mates the energy severely at smaller Ji = J2. We will 
return to this in subsection IIII A 21 

Using DMRG ground state energies, we follow the first 
order transition in the Ji^2-pla'ne. This is shown in the 
inset of Fig. [3] (b). Apart from a very small curvature 
in the immediate vicinity of the transition point on the 
diagonal Ji = J2, the plaquette phase border is composed 
of almost straight lines: J2(i){Ji{2)) ~ 1 for 1 < Ji(2) 
1.5. For values of ^1(2) ^ 1-5, the error on the detection 
of the kink from our numerical data is too large to make 
definite conclusions. While this is identical to previous 
findings in Ref. [l^, our evaluation of the static structure 
factor as in subsection IIII A 3] shows that the first order 
transition is very likely to extend at least up to Ji(2) ~ 5. 

In summary, ground state energy calculations seem 
consistent with a plaquette phase extending throughout 
two strips of width of order unity parallel to each of 
the Ji(-2)-axis, at least up to intermediate Ji{2)- Finally, 
there are no signatures of additional first order transi- 
tions, separating a putative incommensurate and (tt, 0)- 
phase. In view of this 'missing' second incommensurate- 
to-commensurate transition, we will consider also real 
space correlation functions and static structure factors 
m section UnAl 



2. SBMFT phase diagram 

Next we turn to the phase diagram as obtained from 
SBMFT. We use an SU{2) invariant decoupling scheme 
described in appendix I A 21 focusing on solutions with ho- 
mogeneous mean fields. Apart from one Lagrange mul- 
tiplier to fix the local spin, this leads to six bond pa- 
rameters i3„=o,i,2 and ^n=o,i,2j one -B„ and one An for 
each of the three non-equivalent exchange links in Fig. [T] 
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Figure 5: (color online) Non-zero bond mean field parameters 
within the phases of Fig. |4] versus J2 for Ji — 0.01. 



Fig.Hb). 

To conclude, also the SBMFT phase diagram is consis- 
tent with a gapped plaquette phase extending throughout 
two strips of width of order unity parallel to each of the 
Ji(2)-axis, and at least up to intermediate Ji{2)- Within 
this plaquette phase (tt, 7r)((0, 7r))-correlations increase, 
as Ji{J2) increase. Moreover SBMFT shows a (7r,0)- 
phase, similar to the classical case, however with a spin 
gap and without long range order. We return to the 
latter in subsection HiTbI Finally, SBMFT shows no m- 
commensurate phase. 



3. Correlation functions and static structure factor 



Bn refer to triplet, and An to singlet spin correlations. 
Solving the self consistency Eqs. (|A6IA7p cither in the 
continuum limit, or, equivalently minimizing the energy 
of Hamiltonian (|A5[) on sufficiently large finite FFSTs, 
we find the quantum phase diagram shown in Fig. |4] for 
s$ Jl,2 =^ 3. 

First, we emphasize, that the SBMFT solutions in all 
of the parameter space investigated remains gapped. I.e., 
there is no condensation of Schwinger bosons, and corre- 
spondingly no long-range magnetic order (LRO). This is 
to be expected in ID. The 'pitch' vector labels in Fig. [4] 
refer to short range correlations as depicted in Fig. [5l 
which shows a vertical cut through the phase diagram 
of Fig. H close to Ji = 0. In the 'red' phase the AFM 
bond mean fields along the plaquette rungs and the di- 
agonal J2-links are finite, while there are ferromagnetic 
correlations along the Ji-links. In this sense this is a 
(0, 7r)-phase, similar to Fig. [2l The same notion applies 
to the (tt, 0)- and (tt, 7r)-phase. All transitions between 
red, green, and blue phases in Fig. |4] are of first order. 

Fig. [5] clearly shows, that upon lowering J2 the FFST 
continuously evolves into a weakly coupled plaquette 
regime in the red phase. I.e., for J2 ^ 1, the singlet am- 
plitudes Aq oil the plaquette rungs increase up to their 
maximum possible value of 1/2 at J2 ~ 0.25, while the 
inter-plaquette coupling amplitudes jointly decrease to 
zero. Qualitatively similar behavior applies to Ji values 
other than that chosen in Fig. [5] within the red phase 
and within the green phase by interchanging JiO J2 and 

^2,^1 O 52,^1- 

However, as signaled by the grey phases in Fig. 21 and 
from Fig. [5l the SBMFT overestimates the stability of 
decoupled singlet sub-units within the FFST - such as 
the four-spin-plaquette. These grey phases are artifacts 
of the SBMFT which are reached through second order 
transitions. As Fig.[5]shows, SBMFT allows for small but 
finite parameter ranges with only one non-zero and max- 
imized AFM bond mean field, implying that the FFST 
decomposes into a collection of completely decoupled Jq-, 
Ji-, or J2-chains. I.e. in the grey regions the SBMFT is 
incapable to lower the system energy by quantum fiuc- 
tuations between the latter decoupled chains. This is 
the reason for the poor SBMFT ground state energy in 



In this section we turn to the question of a potentially 
incommensurate phase in the quantum case. To this end 
we first look at static real-space correlation functions 

C(r) = (S(r).S(0)), (2) 

where r is a site on the lattice and (. . .) the ground- 
state expectation value. Due to the SU{2) invariance 
of the model, only the correlation function Cz{r) ~ 
{Sz{r)Sz{0)) needs to be considered, which satisfies 
Cz(r) C(r)/3. We will contrast rcsufis from DMRG 
against those from SBMFT. 

SBMFT results are obtained with periodic boundary 
conditions (PBC). For best convergence, DMRG employs 
open boundary conditions (OBC) along the chain. I.e., 
correlations depend on the reference site. To minimize 
edge effects, we have chosen a reference site = (L/2, y) 
in the middle of any of the y = (1 ... 4) equivalent chains 
of the tube. Fig. [6l panels(b), (c), shows C{x) along one 
of those equivalent chains, say r = {L/2 — 1 4- x, 1) and 
= (L/2,1). 

We have focused on three particular values of Ji,2 
as shown in the schematic phase diagram in panel (a). 
Two of them lie regions where both, the classical and 
the SBMFT suggest strongly commensurate correlations, 
and one is shortly above the first order transition of 
Fig. [3l where the classical state is incommensurate. 

Fig. ini^b) evidences clearly commensurate correlations 
along the tube's legs for the regions of the black and 
red open circles in Fig. [G^a) and obviously a remarkably 
good agreement between DMRG and SBMFT. Small de- 
viations between DMRG and SBMFT at the ends of the 
chain are to be expected from the difference in bound- 
ary conditions. We have checked, that the wave vector 
of the commensuration is (tt, 7r)((7r, 0)) for the black(red) 
circles of ^a.) by also scanning along other real-space 
directions on the FFST. Clearly C{x) decays as a func- 
tion of X. While the system sizes for the DMRG are 
too small to extract the functional form of this decay, 
C{x) ~ exp(-a;/C) is found in the SBMFT, where ^ is 
a finite correlation length related with the inverse of the 
energy gap. This is consistent with gapped phases and 
no LRO, as has already been alluded to in section [Tll A 21 
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Figure 6: (color online) (a) Open colored circles: Ji,2 values 
chosen for panels (b) , (c) . PQ refers to weakly coupled plaque- 
tte phase, with predominant (tt, it) (green gradient) or (0, tt) 
(red gradient) correlations. Magenta line: first order transi- 
tion evidenced from DMRG in Fig. E^b, inset) and SBMFT 
in Fig. H (b) and (c) C(x) = {S(a;, 1) • S(L/2, 1)) vs x. Curve 
colors correspond to choices in (a) . Crosses: DMRG for L=30, 
m=300, and OBC. Circles: SBMFT for L=400 and PBC. In- 
set in (b): SBMFT at large distance. 



The situation changes drastically at the blue open cir- 
cle in fig. El^a). Here, DMRG evidences a strongly 
decaying, incommensurate x-depcndence in Fig. [6ljc), 
while SBMFT continues to display commensurate (tt, 0)- 
corrclations, as to be expected from the phase diagram. 
Fig. H This proves, that SBMFT fails to produce the 
proper spin-correlations shortly above the first order 
transition out of the plaquette phase and suggests the 
presence of an incommensurate region also in the quan- 
tum version of the FFST. 

To further corroborate this, we now calculate the static 
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Figure 7: (color online) (p) Letters on yellow background: Ji,2 
choices for panels (a)-(f). PQ refers to weakly coupled plaque- 
tte phase, with predominant (tt, tt) (green gradient) or (0, tt) 
(red gradient) correlations. Magenta line: first order tran- 
sition evidenced from DMRG, Fig. E^b, inset) and SBMFT 
Fig. H Panels (a)-(f): structure factor S(Q) from DMRG 
(L=20, m=300, OBC) for Ji,2 as in (p). Blue, green, red, 
and magenta lines refer to Qy — (0, 1, 2, 3)7r/2. 



structure factor 

^(Q) = ijE^^'^""(S(r)-S(0)), 



(3) 



versus wave vector Q = {Qx,Qy) from our DMRG data 
where = (L/2,1). First we consider a coarse grained 
set of Ji,2. The results are shown in Fig. [71 As labeled in 
panel (a) four values of Ji,2 are taken from regions where 
commensurate correlations are to be expected and two 
out of the vicinity of the first order transition as observed 
in DMRG, Fig.Ei;b, inset) and SBMFT Fig.H Since the 
transverse momentum space of the tube is confined to 
Qy = (0, 1, 2, 3)7r/2 there are four S{Qx,Qy)-^incs for 
each value of Ji,2- 

Fig. El^a) exhibits a flat structure for all Qy modes 
c, which reflects the decoupling of the plaque- 



vs 



y,max 



TT con- 



ttes. Moreover S'(Q) is maximum at Q 
sistent with the singlet ground state on the decoupled 
plaquettes. Figs. [71[d,e,f) show maxima in S{Q) at 
Qmax = ('"',0), (0, tt), and (tt, tt) respectively. This is 
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Figure 8: (color online) S(Q) from DMRG, L=30, m=300, 
and OBC versus Ji = J2 at Q = (i", 0) (blue circles), 
(tt, tt) (red circles) and for the incommensurate Qmax ~ 
±(37r/4, 7r/2) (green circles). Inset: S{n,TT) versus Ji at 
J2 =0. 



consistent with SBMFT in Fig. 3] and also with the classi- 
cal phase diagram in Fig.Hl The small oscillations around 
the maxima are finite size effects. On the finite system 
used for the DMRG calculations, the amplitude of the 
structure factor remains finite at Clmax- From the analy- 
sis up to now, we expect no LRO on the quantum FFST, 
i.e. a finite value of S{Clrnax) for L — c». A proof of the 
latter would require finite size scaling analysis, which is 
beyond our computational reach. 

Figs. [7Ub,c) describe Ji_2 values shortly below and 
above the first order transition of Fig. [3ljb, inset), along 
the line of maximum frustration. Panel (b) contains a 
small modulation in all modes, although the plaquette 
phase is still evident from Qy^max = tt. Panel (c) how- 
ever shows two-symmetric maxima at incommensurate 
vectors with Q^ax = (37r/4, 7r/2), (5^/4,3^/2). While 
the j/-component of these pitch vectors are set by the 
transverse quantization of the momentum space, the x- 
components are set by the quantum correlations in the 
FFST. Very remarkably, these x-components are, up to 
our numerical precision (lO"**), identical to the corre- 
sponding classical pitch- vectors, listed in the enumera- 
tion point [4] in section |lll 

Next, we discuss the DMRG structure factor in a finer 
grained analysis of the Ji^2 plane, along the diagonal line 
of maximum frustration. We have observed the occur- 
rence of Qy^max = TT with flat Qx dependence, charac- 
teristic of the plaquette phase, for Ji = J2 < 0.7 and 
Qmaa; = (j^ i T^)) : Still with vcry flat Qx dependence, for 
0.8 < Jl = J2 < 1 (cf. Figs. ^a,h)). We also flnd 
Qrnax = (tt, 0), signaling a commensurate classical-like 
(tt, 0) phase for Ji = J2 > 1.3 (cf. Fig. [3Ud)). An incom- 
mensurate phase is observed for 1 < Ji = J2 < 1-3, with 
Qmax = (37r/4, V2), (5V4,37r/2). 

In order to describe the extent of such incommensu- 



rate region we show in Fig. [8] S'(Q) for representative 
momenta Q — {tt,tt), (37r/4, 7r/2), (tt, 0) in the range 
Jl = J2 e [0,1.5]. Clearly, ^(Tr, tt) is maximum and 
shows only a small variation in the range Ji = J2 G [0, 1] 
(plaquette phase). At Ji = J2 = 1, the structure factor 
is discontinuous. Following that, and in a small window 
of 1 < Jl = J2 < 1.2, S'(Q) is maximum at the incom- 
mensurate wave vector. In the vicinity of Ji = J2 ~ 1.25 
there is a crossover from incommensurate to commensu- 
rate (tt, 0) correlations. These results can be interpreted 
in terms of a small window of an incommensurate phase 
with a weak flrst order transition into the (tt, 0)-phase 
and a kink in the energy which is too small to be de- 
tected from the DMRG calculations in Fig. [3l 

For reference the inset in Fig. [8] reports S'(7r,7r) along 
the Ji-axis, i.e. J2 = 0, where the structure factor is 
maximum for any Ji > 0. This plot shows a continuous 
increase and no signs of phase transitions in this part 
of parameters space. An identical observation applies to 
5(0, tt) along the J2-axis, i.e. Ji = 0, for all J2 > 0. This 
is consistent with the plaquette phase being adiabatically 
connected with the limit of decoupled chains. 

While the discussion in Fig. [8] is confined to the line of 
maximum frustration, we have performed similar analysis 
along additional lines in the Ji^2 plane. These agree with 
a plaquette phase in strips of width one, both, along the 
Jl-, and J2-axis, as in Fig. [H up to values of J1.2 ~ 5. 
This extends the range obtained from the kink in the 
ground state energy in Fig. [3] and Ref. [H. Moreover, 
incommensurate correlations are observed beyond these 
strips, with Qmax slightly renormalized by quantum fiuc- 
tuations with respect to the classical spiral pitch- vectors 
in the enumeration point [4] in section |TT1 Unfortunately, 
the width of the incommensurate region decreases rapidly 
off from the line of maximum frustration and cannot be 
determined accurately enough. 

To summarize, static structure factor calculations sug- 
gest that the at least close to line of maximum frus- 
tration, the plaquette phase-strips undergo a first order 
transition into an incommensurate phase, the extent of 
which is strongly decreased by quantum fluctuations with 
respect to the classical spiral phase. The transition be- 
tween the incommensurate and the (tt, 0) phase appears 
to be very weakly first order. 



B. Strong leg coupling 

In this subsection we consider the limit of Ji^2 3> Jo by 
explicitly setting Jq = 0. Naturally in this limit, a differ- 
ent normalization is required for the exchange coupling 
constants. We set Ji = 1. For Jo = the FFST is unfrus- 
trated, equivalent to an anisotropic twisted square lattice 
on a torus. We will use quantum Monte- Carlo, based on 
the stochastic series expansion, to study its properties. 
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Figure 9: (color online) a) Finite size scaling for 64 ^ L ^ 
512 of the spin gap for different values of J2. b) L = 00 
extrapolated spin gap vs. J2, as well as spin gap of a spin-2 
chain. Inset: susceptibility vs. temperature for L = 512 and 
3-2. = J\. Symbols: QMC data, solid curve: Fade fit, see text. 
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Figure 10: (color online) Contour plots of the dynamic struc- 
ture factor 5'(Q, to) from QMC & MaxEnt for systems of 
L — 32 vs. uj and Qx at Qy = 0. Panels (a-c): for J2/J1 — 1, 
0.5, and 0.1 at T = Q.25J2. Fanel a) inset: (7r,0)-cut of 
^(Q,^;) with labels referring to panels a)-c). Fanels (d-f): 
for T/Ji = 0.5, 0.25, and 0.1 at J2/J1 = 1. Panel d) inset: 
(7r,0)-cut of S(Q,ijj) with labels referring to panels d)-f). 



The real space arrangement of spins in the classical 
(tt, 0) phase at Ji ~ J2 ^ Jo is that of a spin-2 AFM 
chain. While in the quantum model the total spin per 
plaquette is not conserved, it is nevertheless tempting 
to speculate on a gap similar to that of an actual spin- 
2 AFM quantum chain at J2 = Ji. Additionally, upon 
reducing J2/ Ji^O, the limit of four decoupled chains is 
reached, which is a LLQ and therefore shows no spin gap. 

To test these assumptions, we evaluate the uniform 
spin susceptibility x(^) versus temperature T on sys- 
tems of up to L = 512 plaquettes for J2 = 1, 0.75, 0.5, 
and O.25J1. The case of J2 = Ji is shown in the in- 
set of Fig. IHJb). Obviously, the system has gap. The 
value of the gap is extracted from x(T) by fitting the 
low-temperature behavior for 0.0055 ^ T ^ 0.2Ji to 
X(T) « e-^/^Pl{T)/T, where P^(T) is a Fade approxi- 
mant of order [k,l]. The errors of such fits - for a particu- 
lar choice of the fitted temperature interval - can be made 
less than the QMC's error bars, which are not shown in 
Fig. [9l and arc of order of 10~^. Fig. [9l^a) details the 
finite size scaling of the spin gap for 64 ^ L ^ 512. The 
small oscillations of the data in this plot should not be 
confused with QMC errors or deviations from simple scal- 
ing. Rather they are due to the particular choice of the 
temperature interval for the Fade fit. As is obvious from 
this figure, these oscillations are less than the actual fi- 
nite size corrections. Finally, the main panel of fig. [D^b) 



proves our speculation, namely, the spin gap at J2 — Ji 
is close to that of a spin-2 chainSl and the gap decreases 
monotonously as J2/J1— ^O, where, corresponding to the 
LLQ, A(J2/ Ji = 0) = 0. 

2. Dynamic structure factor 

Continuing on the analogy of a crossover from a gapped 
Haldane-like spin-2 AFM chain to a LLQ for J2/J1 rang- 
ing from 1 to 0, the dynamical structure factor of the 
FFST should show signatures of deconfinement from 
gapped 'magnon'-like modes at J2/J1 = 1 to a two- 
spinon continuum as J2/J1— >0. 

To analyze this, we investigate the dynamic struc- 
ture factor S'(Q,w) at frequency w, which we obtain 
from MaxEnt analytic continuation of imaginary time 
dynamic structure factor 

^(Q'^) = izE^''^""(S(r,r)-S(O,0)), 

r 

evaluated by QMC (see appendix IA3p . This is shown 
in Fig. [TUl In all of these plots Qy = 0. The absolute 
scales on all panels of this figure are adjusted to ensure 
approximately identical extent of the spectra along the 
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y-axes, which allows to compare the width of the spectral 
contours. Turning to Fig. [TCT a-c). we first note that all 
three contour plots display a certain broadening due to 
the finite temperature T = 0.25Ji. We return to this 
in Fig. fTOFd-f). Apart from that, at Ji = J2 the figures 
show a rather sharp magnon-like mode, similar to the 
spectra of integer-spin Haldane chainS j^^i^^ accompanied 
by a marked loss of spectral weight as Qx — 0, which is 
also a typical feature of integer spin chains <^ As J2 — !• 0, 
the spectrum starts to broaden in the vicinity of = tt, 
resembling a shape very similar to that of the spinon 
continuum of the spin-1/2 AFM Heisenberg chain^i^ - 
exactly as anticipated. The inset in Fig. [TOFa) details, 
that although the finite temperature maximum of the 
dynamic structure factor does not have to coincide with 
the spin gap, it nevertheless decreases similar to the latter 
with respect to J2I Ji- 

Figs. [TUTd-f) list the temperature dependence of 
^(Q,^) for J2 = Ji- First, these panels clarify, that 
T = O.25J1 is a reasonable compromise between finite 
size effects at L = 32 and thermal broadening, i.e., for 
T = O.lJi the line broadening is already less than the 
finite-size level-spacing. Furthermore, the inset llOf d) col- 
lects cuts at Qx = TT, which demonstrate a rather strong 
temperature dependence of the zone-boundary modes of 
the FFST for J2 ~ Ji > Jo- This might be of interest in 
the context of similar observations^^ for four-spin tube 
compound Cu2Cl4-D8C4S02. 

In conclusion, and even though the on-plaquette to- 
tal spin is not strictly conserved in the quantum case, 
the FFST shows some features remarkably similar to an 
AFM chain with 5 w 2 at Ji = J2 ^ Jo, as well as a 
magnon-spinon deconfinement as J2/J1 0. 



IV. DISCUSSION AND CONCLUSIONS 

To summarize, we have studied the quantum phases 
of a frustrated spin-1/2 four-leg tube using a variety of 
techniques: density matrix renormalization group. Quan- 
tum Monte-Carlo, Schwinger boson mean field theory, 
exact diagonalization and series expansions. Our main 
results arc outlined in the tentative quantum phase dia- 
gram Fig. [TTJ This figure should be contrasted against 
the tube's phase diagram in the classical limit, i.e. Fig.[2l 
While all phases in the latter are long range ordered, none 
of the quantum phases are. 

The point Ji^2 = hosts a gapped system of decoupled 
plaquettes, while at the asymptotic points Ji(2) ^ 00, 
■^2(1) = the spin tube degenerates into decoupled spin- 
1 /2 chains in a Luttinger liquid state. The phase diagram 
is symmetric with respect to interchanging Ji -;-> J2. On 
either of the two axes Ji{2) = 0, the system is unfrus- 
trated, the inter- leg coupling is relevant, and a spin gap 
opens. This unfrustrated weakly coupled chain regime 
is known to be adiabatically connected to that of the 
weakly coupled plaquettes. 

Turning on the frustrating exchange, our results arc 



OO" 



A 

o 

_l" 
I 

5 



-deconf. -AFM "spin 

N) 
O 



J2 



2 + 
1 




gapped 

(^,0) 



^-1"'ord 
I ^weak 1^* ord 




o 
o 



PQ 
— I— 



H 1- 



1 2 



Ji 



4-LLQ 

' 00 



Figure 11: (color online) Quantum phase diagram of the 
FFST 



consistent with the weakly coupled plaquette regime to 
survive along two strips (red and green in Fig. [TT|) of 
width of order unity, parallel to each of the Ji(2)-axis, 
at least up to Ji(2) ~ 5 Jq. The system remains gapped 
in this region. Accordingly, our analysis of correlation 
functions exhibits exponential real space decay. Consis- 
tent with series expansions around Ji_2 — 0, the static 
structure factor obtained from density matrix renormal- 
ization group evolves smoothly from a fiat plaquette sig- 
nature around PQ in Fig. [11] into a peaked commensu- 
rate behavior along the red/green strips, parallel to each 
axis. The peak locations are consistent with short-range 
correlation remnants of the long-range order present in 
the classical limit of the tube this region. As for the 
unfrustrated four-leg tube, we expect no quantum phase 
transition while increasing Ji(2) — ^ 00 parallel to the axis 
within these strips until the Luttinger liquid fixed point 
is reached (zig-zag marks in Fig. [TT|l . 

Perpendicular to the Ji,2-axis the plaquette regime 
is terminated by a line of first order transitions evi- 
denced by those of our techniques able to detect ground 
state energy level crossings. The critical lines emerge 
approximately from the point of maximum frustration 
•-'0.1,2 = 1 and run parallel to the Ji 2 axes (magenta line 
in Fig. Ill|) . The numerical precision, locating the level 
crossing along the borders of the PQ strip, decreases away 
from Jo. 1.2 = 1, indicated by the doting of the magenta 
line. 

Beyond the first order critical line, close to the point 
of maximum frustration, Jo. 1,2 = 1, DMRG shows that 
the plaquette phase turns into a gapped phase with short 
range incommensurate correlations (IC, blue in Fig. Ilip. 
analogous to the spiral phase which is found in the clas- 
sical limit of the tube in this regime. Along the diagonal 
1 < Ji = J2 < 1.3, the static structure factor shows 
a maximum approximately at the pitch vectors of the 
classical spiral phase. Off the diagonal, the maximum 
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of the static structure factor is slightly shifted from the 
classical values. Increasing the inter-plaquette coupling, 
around the line Ji ^ J2, the incommensurate quantum 
phase terminates with a very weak first order transition 
into a gapped commensurate (tt, 0) phase, labeled by the 
thin black line in Fig. [TT] In contrast to the PQ (tt, tt) 
and (0,7r) region, the overall extent of the incommensu- 
rate region in the quantum case is strongly reduced as 
compared to that of the classical spiral phase. 

For Ji_2 ^ Jo, the system can be considered as approx- 
imately unfrustrated. We have investigated this regime 
by Quantum Monte-Carlo along the line < J2/J1 < 1, 
setting Jo = (black line emerging perpendicularly from 
the point LLQ in Fig. [TT|) . Here calculations of the uni- 
form susceptibility show the tube to have a gap very 
close to that of AFM spin-2 chains at Ji = J2, while 
for J2/J1 the gap decreases to zero as expected 
for approaching the Luttingcr liquid state. Evaluat- 
ing the dynamic structure factor, and consistent with a 
crossover from a 'Haldane-like AFM spin-2 chain' behav- 
ior at Ji ~ J2 to a LLQ at J2 = 0, we observe a decon- 
fincment of the excitations turning from sharp 'magnon' 
modes into a spinon continuum as J2/J1 ^ 0. 

Finally, and due to numerical limitations in our study, 
it remains an open issue if the quantum IC and PQ 
regimes extend beyond Ji(2) ~ 5 at J2(i) ~ 1. In this 
context wc cannot conclude from our study whether the 
PQ and (tt, 0) phase remain adiabatically disconnected 
in the quantum case or not. 
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Appendix A: Techniques 

For completeness, this appendix provides some details 
and references to the methods we use in this work. 



1. Series expansion 

Our SE calculations start from the limit of isolated 
plaquettcs. To this end we decompose the Hamiltonian 



of the FFST into 

H = Ho + V{JuJ2). (Al) 

where Hq represents decoupled plaquettes and V^( Ji, J2) 
is the part of Hamiltonian that connects plaquettes via 
Ji, J2 couplings. 

It is simple to show, that each plaquette has four 
equally spaced energy levels, which in turn renders the 
levels structure of Jfp to be equidistant. This allows to 
sort the spectrum of Hq in a block-diagonal form, where 
each block is labeled by an energy quantum- number Q. 
In this way, Q=0 represents the ground state {vacuum), 
i.e., all plaquettes are in the state of minimum energy. 
Q=:l sector is composed by states obtained by creating 
(from vacuum state) one-elementary excitation {particle) 
on a given plaquette, and so on. It is clear that Q > 2 
will be of multiparticle nature. 

In general the action of V^(Ji, J2) mixes different Q- 
sectors, so that the block-diagonal form of Jfg is not 
conserved in H . However, it has been showr^ that for 
the present type of Hamiltonians it is possible to restore 
block-diagonal form by the application of continuous uni- 
tary transformations, using the flow equation method of 
Wegner— . It basically consists in transforming H onto 
an effective Hamiltonian iJoff which is block-diagonal in 
the quantum number Q. This transformation can be 
achieved exactly in terms of a SE in J1.2 leading to 

^off ^ Hq + ^ J" "^J^Cn,m ■ (A2) 
n,0<rn< n 

Here C„^„i are weighted products of terms in V^(Ji, J2) 
which conserve the Q-number, with weights determined 
by recursive differential equations (see Refi^ for details) . 

Due to Q-number conservation several observables can 
be calculated directly from HcS in terms of a SE in Ji_2- 
For systems with coupled spin-plaquettes continuous uni- 
tary transformations SE has been used for one^^, two^^Ti^^ 
and three^ dimensions. For the present model we have 
performed 0{7) and 0{6) SE in J1.2 for ground state en- 
ergy {Q — 0) and for Q = 1,2 sectors, respectively We 
refer for technical details about the calculation to Ref J^. 



2. Schwinger bosons 

Schwinger bosons^ are used to represent spins at site 
I via spinfuU bosons b'j^J , with a =ti or ±1, through 

•5;" = ^T,,,i.bl^<^'ilAu^ where are the Pauli ma- 
trices and a = x,y,z. The Hilbert space dimension 
of spin-5 multiplets is enforced through the constraint 
^la^ia~'^^ ■ tcrms of Schwiuger bosons, the ex- 
change interaction can be written as^^*^ 

Si ■ S,n = : bIjiBi„^ : -^!mAm ' (^3) 

with the bond operators Bj^ = ^J2o- ^L^mo- ^^'^ ^'"i 
= crbiabm-(T and normal ordering ::. Eqn. (|A3P 
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has been used for various SU{2) invariant and large N 
factorization schemes^"—. We foUo w^^'^^ and introduce 
the bond mean fields Bim, = {Bim) and Aim = {Aim), 
accounting for ferromagnetic (FM) and AFM correlations 
on equal footing. For the FFST we focus on homogeneous 
mean fields, implying six parameters 



Br, 



(A4) 



where n = 0, 1, 2 corresponds to the three exchange links 
1"; — r,„=r„=Ry, R-c, R^; + Ry. Fourier transformation, 
bia = X^k '^''""'^kcr/v^, leads to a bilinear mean field 
Hamiltonian, which can be diagonalized by standard Bo- 



goliubov transformation, i.e. h^^ 
with — = 1 yielding 



MFT 



kcr 



= E^k(aL«k.+i)+y:^4iA/ 



2N\{S- 



2^' 



(A5) 



where = [7b (k)^ 



7A 



(k)^]^/^ is the quasiparticle 



dispersion with 7a (k) = JnAn sin(k • r„) and 7s (k) 
= ^„ JnBn cos(k- r„) — A. We assume Bn, An to be real. 
A is a Lagrange parameter to enforce the constraint on 
the average. Selfconsistency, i.e. d{HMFT) /dx^O, with 
x=An, Bn, and A leads to 



A[Bl 



-E 



7A[i3](k) sin[cos](k • r„) 



2A^ 



1 >p 7B(k) 
2A^ ^ Ek 



(A6) 
(A7) 



where eqn. (|A6p yields six equations for An and Bn, by 
replacing terms with their square bracketed successors. 

To obtain A„,_B„, and A we use two numerical ap- 
proaches: (i) we solve eqn. (|A6IA7p in the thermody- 
namic limit, and (ii) we minimize the vacuum energy of 
eqn. (jASP with respect to An , Bn , and A on large finite 
lattices with N ^ 10^ sites and periodic boundary con- 
ditions. The results from both approaches agree. 

In the present work we set = 1/2 and study the 
ground state energy, the quantum phases, and the spin 
correlation functions arising from An,Bn, and A. 



where Hi , 



1/2. 



S,,,.^St,,.)/2 are spin diagonal and off-diagonal bond op- 



S^) refers to the 



Kb) 

orators between sites i, J. \a) = \Sf 

basis and Sn = [ai, &i] [02, ^2] ■ • • ^n] is an index for 
the operator string Ofci -^a.k,bk- This string is Metropo- 
lis sampled using diagonal updates which change the 
number of diagonal operators Hi^bk in the operator string 
and directed loop updates which perform changes of the 
type Hi b^. ^^2,&fc- For unjrustrated spin-systems the 
latter update comprises an even number of off-diagonal 
operators -^2,6^ , ensuring positivencss of the transition 
probabilities. 

Evaluation of the transverse dynamic structure factor 
with QMC is performed in real space i,j and at imagi- 
nary time r following Ref4ii 



St{p)Sj{p + m)\ , (A9) 
/ w 

where {. . )w refers to the Metropolis weight of an op- 
erator string of length n generated by the stochas- 
tic series expansion of the partition functio n^^'^^ , and 
p, m are positions in this string. Analytic contin- 
uation to real frequencies follows from the inversion 
of S^{c\,t) = i duS^{(\,ijj)K{ijj,T), with a kernel 
A'(tj,r) = e-™ + e-(/3-r)cu p ^ -^^j,^ S'±(q,T) = 
Eae*'*""^<^.o(r)/iV. 

The preceding inversion is performed using the max- 
imum entropy method (MaxEnt), minimizing the func- 



tional Q = X /2 — ao 



43,44 



Here x refers to the co- 



variance of the QMC data to the MaxEnt trial-spectrum 
Sa±i<i,i^)- Overfitting is prevented by the entropy a = 
"^ailq, w)ln[5Qj_(q, w)/m(w)]. We have used a fiat 
default model m{u}), matching the zeroth moment of the 
trial spectrum. The optimal spectrum follows from the 
weighted average of Sa±{q,uj) with the probability dis- 
tribution P[a|S'(q, r)] adopted from Ref4^. 



3. Quantum Monte-Carlo 

We employ the stochastic series expansion (SSE)^"— , 
which is based on importance sampling of the high tem- 
perature series expansion of the partition function 

a n S„ ' k=l 



4. Exact diagonalization and density matrix 
renormalization group 

All DMRG and ED calculations employ the open 
source packages ALPSi^ and SPINPACK^. We refer to 
their documentation. In DMRG specifications, m refers 
to the number of states kept during sweeps. 
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